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Abstract 

We consider the fist order, gradient-flow, description of the scalar fields coupled to 
spherically symmetric, asymptotically flat black holes in extended supergravities. 
Using the identification of the fake superpotential with Hamilton's characteristic 
function we clarify some of its general properties, showing in particular (besides re- 
viewing the issue of its duality invariance) that W has the properties of a Liapunov's 
function, which implies that its extrema (associated with the horizon of extremal 
black holes) are asymptotically stable equilibrium points of the corresponding first 
order dynamical system (in the sense of Liapunov). Moreover, we show that the fake 
superpotential W has, along the entire radial flow, the same flat directions which 
exist at the attractor point. This allows to study properties of the ADM mass also 
for small black holes where in fact W has no critical points at finite distance in mod- 
uli space. In particular the W function for small non-BPS black holes can always 
be computed analytically, unlike for the large black-hole case. 
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1. Introduction 



It is well known that in A^-extended supergravity based on symmetric coset manifolds G/H the 
dynamics of extremal (spherically symmetric, asymptotically flat) black holes is encoded in a "fake" 
superpotential function W [IJ [2| [3l 131 [6^1 which. for large black holes, is entirely specified by the duality 
orbit [8] of the dyonic charge vector V — (p A , q\) (A = 1, • • • , ny) and the asymptotic values at radial 
infinity of the scalars of the theory: (<%) £ G/H. 

It has not been appreciated enough that many properties of the W function are not only true at the 
horizon of regular extremal black holes, defined by critical points 0* of W: 



dW 



= 0, (1) 



where space-time is AdS2 X S2, but in fact they are valid on the entire radial flow and in particular at 
spatial infinity where space-time is flat. All these properties naturally follow from the identification of the 
M^-function with Hamilton's characteristic function [3J[H] of an autonomous Hamiltonian system|f|. In [5], 
the fact that the radial evolution of scalar fields and the metric of spherically symmetric, asymptotically 
flat solutions, is described by an autonomous Hamiltonian system, was used to show that the problem 
of defining a first order description of black holes in terms of a superpotential W is equivalent to a 
Hamilton- Jacobi problem: The W function is a solution to the Hamilton- Jacobi equation associated with 
the Hamiltonian system and defines a system of first order gradient-flow equations for the scalar fields. 
From the identification of the superpotential W with Hamilton's characteristic function, some important 
general properties follow: 

• W is a positive definite function on the moduli space; 

• The derivative of W along the flow, moving from the horizon to radial infinity, is always positive; 

• W for extremal solutions is duality invariant. 

The latter property was originally conjectured in [2] and later proven in [3J. Eventually, in [H [SJ [B] the 
explicit construction of W in terms of duality invariants was completed. 

As we shall prove, the above properties of W, in the presence of a critical point 0* of the first order 
system, promote W to a Liapunov's function (see for instance allowing to make precise statements 

about the asymptotic stability of </>* (namely that <fi* is not just an attractive equilibrium point, but 
also stable), with no need of computing the Hessian of the potential. The existence of W, even in a 
neighborhood of the critical point, provides an alternative (and more powerful) characterization of its 
attractiveness and stability properties. 

Of particular interest are orbits of extremal, large, black holes[f| in which the critical points are not 
isolated. This feature is related, in the symmetric models, to the existence of flat directions of the scalar 
potential V [Hj. We shall prove in full generality, using the duality invariance of W, that W and V 
have the same symmetry properties, and thus that they also have the same flat directions. These flat 



1 The idea of "fake" supersymmetry was first introduced in [7]. 
2 See also jTU] for related independent works. 

3 Large black holes are solutions for which a certain (quartic) duality-invariant expression of the charge 
vector V, called A('P) does not vanish. For small black holes, on the other hand, Ii{V) = 0. A definition 
of Ii(V) ,and its G-invariant form for symmetric geometries G/H is summarized in subsection 2.1 
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directions, which have an intrinsic group theoretical characterization in terms of the stabilizer of the 
duality orbit of the quantized charges [8], are also a feature of the central and matter charges. 

An other consequence of the general properties of W is that the functional form of W{i n ,I^, where 
i n are the H- invariant combinations of the moduli <ff and charges V, can also be calculated for I4 = 0, in 
which case the classical horizon area vanishes and eq. |T]) has no solutions (in the interior of the moduli 
space). More precisely, for I4 = 0, eq. (JJ) has a runaway solution W — at the boundary of the moduli 
space where some <ff — > 00 [13] , 

It is the aim of this note to further specify general properties of the W function for large and small 
black holes, such as their moduli spaces and symmetries. Moreover, depending on the number N of 
supersymmetries, W(Ii = 0) can be obtained by a suitable limit of large black hole solutions (where 
I a 7^ 0), in such a way that W is always given by a calculable algebraic function of the H- invariants. The 
way the limit is performed also allows us to simply understand the interplay of BPS properties of small 
black holes versus large solutions. 

As a byproduct, our analysis shows that I4 — black holes in J\f — 8 are always BPS (distinct in 
three orbits with different fractions of supersymmetry), while in M — 4,2 theories I4 — black holes 
can be non-BPS. Nevertheless their W function in this case can also be analytically computed from a 
non-BPS large black hole with vanishing central charge (Zab = 0, I4 > 0). 

The paper is organized as follows: In section 2 we review and outline the above mentioned properties 
of the W supcrpotential. In particular wc give the general form of W for extremal solutions and address 
the issue of asymptotic stability of the critical points on W, by showing that W can be identified with a 
Liapunov's function. We also prove that the potential V, the superpotential W, together with the central 
and matter charges, have the same flat directions defined by the G-orbit of the quantized charges V . In 
sections 3, 4, 5 we analyze small black holes for M — 8, 4 and 2 respectively. Three appendices containing 
derivations of properties of the W-function discussed in the text, together with other mathematical 
details, are included. 



2. Some General Properties of the W Function 



Let us review some general properties of the fake superpotential W associated with [/-duality orbits 
of static, extremal, asymptotically flat black hole solutions in an extended supergravity theory with a 
symmetric scalar manifold ^# sca z 



H ■ 



Let us consider an extended supergravity describing n real scalar fields, spanning the manifold ^# sca i 
and ny vector fields A^. The ansatz for the metric and the vector field strengths F^ v , for the kind of 
black holes we are considering, is: 



ds 2 



J1U 



dt 2 



-2U 



dr 2 1 . ~ 

— +-2 v<p 



F = 



dx^ A dx u 



„2U„ 



•Ml 



■sm(9)dip 2 ) 
")-VdtAdT + T sin(9)d0 A dip, 



(2) 



where the coordinate r = — 1/r runs from 0, at radial infinity, to —00 at the horizon, where e u ^ vanishes. 
The scalar fields are taken to be functions of r only: (ff — 4> t (t). The magnetic field strength Ga^ in (0) 



is defined, as usual, as: G 



oc e 



8£/5F A a , C begin the Lagrangian of the theory. The last equation 



in ^ is written in a manifestly symplectic covariant form, namely as an equality between two 2 ny 
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dimensional symplectic vectors, where Cmn, M, N = 1, . . . , 2ny is the Sp(2riv r ,K)-invariant matrix: 

The vector V = (p A , q\) consists of the quantized electric and magnetic charges. Finally the 2 ny x 2 ny 
symmetric, negative defined, symplectic matrix M.{<ff)MN = — (LL t )mjv, L(</> r ) being the ^ sca i coset 
representative in the fundamental of Sp(2ny,R), can also be written in the familiar form [16 : 

™ - (™ -r ) • 

where /as = Im(^V)AS < is the vector kinetic matrix while i?A£ = Rc(A/")as defines the generalized 
theta-term. 

Once the electric and magnetic charges of the solution are assigned, the radial evolution of the n + 1 
fields U(t), r (r) is described by the effective action: 

Seff = J £ ef fdT = J (u 2 + ^G r M<P r <P s + e 2U V(0,V)^j dT, (5) 

together with the Hamiltonian constraint, representing the extremality condition^: 

H ef f = U 2 + ^G rs (cl ) )^^-e 2U V(4>,V)^0, (6) 

the effective potential being given by V(<p, V) = — ^ "P T .M(0) V > and the dot represents the derivative 
with respect to r. The radial evolution of the n + 1 fields U(t), </>' (t) in the solution admits a first order 
description [TJ[21[3] in terms of a fake superpotential W((j),P): 

U = e u W , cj) r = 2 e u G rs — . (7) 
Y dcf) s w 

If we interpret the fields U(t), <P t (t) as coordinates of a Hamiltonian system in which the radial variable 
plays the role of time, the first order description ([7]) is equivalent to solving the Hamilton- Jacobi problem 
with Hamilton's characteristic function 

W(U,(f)) = 2e u W(<f>). (8) 

Indeed, in terms of W(<j),V), the Hamilton- Jacobi equation has the form: 

dW dW 

which can also be derived from the Hamiltonian constraint d6l) using We are not interested here in 
the most general solution to (j9]), nor to address the issue of its existence (see [3] for a discussion on this 



4 For non extremal black holes the value of the Hamiltonian on a solution coincides with the square of 
the extremality parameter. Notice that the Hamiltonian is not positive definite, being expressed as the 
difference of a "kinetic" and a positive "potential" term (this is in turn due to the fact that the role of 
the time variable is played by a spatial coordinate t). As a consequence of this we can have non-trivial 
solutions on which the Hamiltonian vanishes. These correspond to the extremal black holes. 

5 In the case of non-extremal solutions the Hamilton- Jacobi equation reads (|^) 2 + 2 G rs (4>) |^ = 
4 e 2U V + 4 c 2 , and the corresponding first order equations have the form V — \ f^j-, 4> r — G rs (cf>) If 
c ^ however, as it is apparent from the Hamilton- Jacobi equation, the dynamical system can have no 
equilibrium point = = 0. 
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point). We are interested, instead, in the W functions associated with classes of extremal solutions whose 
general properties are in principle known. They are completely characterized by the set of quantized 
charges V and the values of the fields at radial infinity: 

U(t = 0) = , 0> = O) = 0S. (fO) 

We shall therefore simply denote them by: U = C/(r;0 o ) an d (jf = <fi r (t; tfio). The ADM mass and the 
scalar charges at infinity are given by: 

M ADM {<f> Q ,V) = U( T = 0) = W(<f> ,V), 

dW 

TT{<f> Q ,V) = <t) r(T = {)) = 2G rs {4>v)—{^V). (II) 

Regular (large) extremal black holes have finite horizon area Ah and thus near the horizon (t — » — oo) 
e u has the following behavior: e~ 2U ~ ^-t 2 , where Ah = Ah(V) is a function of the quantized 
charges only. In fact V transforms under duality (see subsection 2.2) in a symplectic representation of 
G and Ah, as a function of V , is expressed in terms of the quartic invariant of G in this representation: 
Ah(V) = 4 7r \J\I&(P)\ (here we use the units c = h = G = 1, so that the Plank length is one.). Using 
eq.s we see that W computed on the solution evolves, in the near horizon limit, towards J 
As far as the scalar fields are concerned, due to the attractor mechanism some of them are fixed at the 
horizon to values which are totally determined in terms of the quantized charges, while other scalar fields, 
which are flat directions of the potential, are not. That is, in the presence of flat directions, in the near 
horizon limit r — > — oo the non flat scalars evolve towards values which are totally fixed in terms of 
quantized charges, while the flat directions still depend, in general, on the boundary values $j taken at 
radial infinity (t = 0). Since only the scalars parametrizing the flat directions may depend at the horizon 
on </>q, the near horizon geometry, which is determined in terms of the potential, will only depend on the 
quantized charges, consistently with the attractor mechanism. Summarizing, for large black holes, we 
have: 



lim e- 2U = VMPT\T 2 , hm f(r) = # , 

T — y — OO T — ^ — OO 



lim W 2 (<I>(t;M,V) = W 2 (<f>*,V) = V(<f>*,V) = ■ 

T— > — OO 

Small black holes are characterized by vanishing horizon area, i.e. by quantized charges for which Ia(V) = 
0. For t — > — oo the warp factor has the following behavior: e~ 2U ~ r a , a < 2. In the same limit scalar 
fields typically flow to values which are at the boundary of the scalar manifold. Either for large or for 
small solutions, from the first of ([7]) we deduce the following boundary condition for W: 

lim e u ( r ' M WMt-^o),?) = lim {7 = 0. (12) 

r— V — oo r— > — oo 

This allows us to write W(4>,P) for the two kinds of solutions in the following form (see [3]): 

r° 

W(<f> ,V) = / & tu ^^V{<f>(T](M,V)dr. (13) 



It should be stressed that the above expression allows to write the W function for a given class of solutions 
as a free function of the point (f>o on the scalar manifold and of the quantized charges: Given a charge 
vector V and a point = (<%) in ^ sca i, the corresponding value of W is given by the integral over r of 
e 2U V, computed along the unique solution originating at infinity in <pQ. 

2.1. I4 Invariant for Af = 2,4,8 Supergravities 
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In Af > 2 theories and in Af = 2 theories based on symmetric spaces for the (vector multiplet) scalar 
fields, the entropy area law reads (the Boltzmann constant ks being one in our units): 

S = j=tvM, (14) 

where, as anticipated in the previous section, IiiV) is a certain quartic invariant of the dyonic charge 
vector V and depends on the particular theory under consideration. Since Ii(fP) is moduli-independent, 
it can be expressed either in terms of the quantized charges V or in terms of the (dressed) central and 
matter charges Zab{4>i "P)i ^i{<j>, V) (see subsection 2.3 for a precise definition of the latter). For our 
convenience we recall here the actual form of I^ifP) in terms of the central and matter charges. 

For Af = 2 theories, based on special geometry, we can define five iJ-invariant quantities i n , as follows 

®: 

ii = ZZ, 



j > 



k = ^Re(ZN 3 (Z,)) , 
U = ~Im(ZN 3 (Z,)) , 

«5 = g u Cijk Cjj^, Z Z ' Z-* Z^ , 

where N 3 (Zi) = djk Z l Z 1 Z k , Zi = DiZ and Z % = g li Zi. In terms of these quantities the quartic 
invariant reads: 

h - (ii-i2) 2 +4,u-i 5 =h(V), (15) 

where, as anticipated in the previous subsection, V transforms in a symplectic representation of G and 
l4{V) is the only non- vanishing invariant quantity built out of the charge vector. Note that, for the 
quadratic series (C^k = 0) we have: I4 = /|, where I 2 = \ii — 

For Af = 4, we can define two SU(4) x SO(n) invariants: 

1 „ —ab 
2 



Si — — Zab Z — Zi Z j S 1 ^ . 



S2 = j e ABCD Zab Zcd — Zj Zj S IJ , 

in terms of the central charges Zab = ^Z^a-, A,B= 1, . . . , 4, and the n matter charges Zj, I = 1, . . . , n. 
Then the unique quartic G — SL(2,R) x SO(6, n)-invariant reads: 

jf= 4) (P) = Sf-lSzW (16) 

and the black hole potential is: 

V^(^T) = \z AB Z AB + Z^ . (17) 

Finally, in the Af — 8 theory the Cartan G = E7(7)-quartic invariant is given by the expression |18j : 
lf^\v) = Tr[(ZZ t ) 2 ]-[Tr(ZZ t )] 2 +8Re[P/(Z)], (18) 



where Z = (Zab) = — Z T , A, B — 1, . . . , 8, is the complex central charge matrix [19]. In terms of the 
four skew-eigenvalues Zi, i = 1, . . . , 4, of Zab, I^ -8 " 1 reads: 

4 

/f=V) = E l z « I* - 2 5^N 2 |%- 1 2 + 4 (Z1Z2Z3Z4 + ziz 2 z 3 z 4 ). (19) 
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The black hole effective potential has the following form: 

V^H^V) = \z AB Z AB = £N 2 ■ ( 20 ) 

In any extended supergravity, BPS solutions are described by W = \zh\, where Zh is the highest skew- 
eigenvalue (i.e. eigenvalue with highest modulus) of the central charge matrix Zab (for N — 2, Zab — 
Z 6ab, A,B = 1, 2, and Zh = Z). Therefore it is also true that: 

V = \z h \ 2 +2G rs d r \z h \d s \z h \. (21) 

If however V is not in a BPS orbit, the flow defined by W — \zh\ does not correspond to a physically 
acceptable solution and a different VF-function should be used. 

In particular, in the N = 8 case for non-BPS configurations the corresponding VU-function satisfies 
the following inequalities: 

\z h \ 2 <W 2 <4\z h \ 2 , (22) 

the lower bound being saturated only for BPS solutions. The upper bound originates from the general 
property: W 2 < V < 4|z/j| 2 . For non-BPS large black holes, it can be proven that, at the attractor 
point, \z%\ = p — \zh\ and the upper bound is saturated: W = 2 p. 

2.2. The W Function and Duality 

It is known that the on-shell global symmetries of an extended supergravity, at the classical level, are 
encoded in the isometry group G of the scalar manifold (if non-empty), whose action on the scalar fields 
is associated with a simultaneous linear symplectic action on the field strengths F A and their duals G\- 
This duality action of G is defined by a symplectic representation D of G: 

( (j) r ->■ <j) r ' = g*(j) r 

where g* denotes the non-linear action of g on the scalar fields and D(g) is the 2n v x 2n v symplectic 
matrix associated with g. The matrix Ai{4>) transforms under G as follows: 

M{g*ct>) = D{ g y T M{4>) Dig)- 1 . (24) 

A duality transformation g 6 G maps a black hole solution U(t), </>' (t) with charges V into a new solution 
U'{t) = U(t), <p' r {T) = g-k<jf{T) with charges V = D{g)V. More specifically, if U(t), (jf{r) is defined 
by the boundary condition 4>$ for the scalar fields, U'(t) — U(t), ,r (r) is the unique solution, within 
our class, with charges V defined by the boundary condition 4>' = g * <p 

(U(t;<P ) fU'(T;g*(f>o) = U(T;<t>o) 

geG : I 0( r; fa) — > <<//(t; g * fa) = g * 0(r; fa) . (25) 
(V [v' = D(g)V 

Using eq.s (fMf and (|2^|) . we see that the effective potential is invariant if we act on (jf and V by means 
of G simultaneously: 

V(<f>,P) = V{g*ct>,D{g)V). (26) 

This implies that V, as a function of the scalar fields and quantized charges, is G-invariant. From this 
property of V it follows that the effective action ([S]) and the extremality constraint ^ are manifestly 
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duality invariant. Let us show now that the W function shares with V the same symmetry property (|26[) . 
namely that it is G-invariant as well: 



— oo 


e" 

-oo 




WfaP) = W(g*<P,D(g)V). (27) 
This is easily shown using the general form (fT3]l and eq.s (|2"5|) : 

W{g*<t> ,D{g)V) = [ e 2U '^°*^V(cf>'(T;g*cb Q ),D(g)V)dT = 

e 2U(r- M y^r-^^dr = W(4> ,V) . (28) 

Being the ADM mass expressed in terms of W, see eq. (fTTj) . it is a G-invariant quantity as well: 

M ADM (<Po,V) = M ADM (g*4>o,D(g)V). (29) 

Extremal black-holes can be grouped into orbits with respect to the duality action (f2"5j) of G. These 
orbits are characterized in terms of G-invariant functions of the scalar fields and the quantized charges, 
which are expressed in terms of iJ-invariant functions of the central and matter charges. One of these is 
the scalar-independent quartic invariant Ii(V) of G which defines the area of the horizon for large black 
holes. Small black holes, on the other hand, belong to the orbits in which l^(P) = 0. 

2.3. The Issue of Stability: Asymptotic Stability of the Critical Points 



Let us notice, from eq. (|T3|, that W is always positive definite, since the effective potential is. 
Moreover its derivative along the solution 4> r (r) is positive definite as well (except in where it vanishes): 

= &d r W = e~ u G rs U)& j) s >0. (30) 

dr 

We see that, if 4>* is isolated, W has the properties of a Liapunov's function and thus, in virtue of 
Liapunov's theorem, 0* is a stable attractor point (we refer the reader to Appendix [Bl for a brief review of 
the notion of asymptotic stability in the sense of Liapunov and of Liapunov's theorem, see also standard 
books like [11]). This conclusion extends to models based on a generic (not necessarily homogeneous) 
scalar manifold: The very existence of the VF-function (i.e. of a solution to the Hamilton- Jacobi equation) 
even just in a neighborhood of an isolated critical point is enough to guarantee asymptotic stability 
of <^», and thus that the horizon is a stable attractor. Let us emphasize that in this case we need not 
evaluate the Hessian of the potential on 0* . In other words the (local) existence of W can be taken as 
an alternative and more powerful characterization of the attractiveness and stability properties of the 
horizon point <f>* . 

There is a class of large extremal solutions, however, in which the critical points, defining the near- 
horizon behavior of the scalar fields, are not isolated but rather span a hypersurface C of the scalar 
manifold. This is the case of the non-BPS solutions with I4 < in the symmetric models. As we are 
going to show below, in full generality, the existence of this locus of critical points is related to the 
existence of nf < n flat directions tp a , a — 1, . . . , nt, of both the scalar potential V and the W function. 
The critical hypersurface C has in this case dimension rif and is spanned by ((p a ). As far as the global 
behavior of the flows is concerned, the analysis of the simple STU model (see [T7] for a discussion on this 
point) suggests a general property: The scalar manifold can be decomposed in hypersurfaces ^#( Q ) of 
dimension n — nf which intersect the hypersurface of critical points C in a single point 4>*\ a characterized 
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by fixed values tp a of the flat directions. The hypersurfaces ^( a ) have the property of being invariant 
with respect to the flow, namely that, choosing the initial point cf>o on a given ^# Q , the entire flow will 
be contained within the same hypersurface. Within each ^#( a ) the critical point 4>*\ a is isolated and 
Liapunov's theorem applies, implying it is asymptotically stable or, equivalently, a stable attractor. 

2.4- The Issue of Flat Directions 

Let us denote by Go C G the little group (or stabilizer) of the orbit of the quantized charges V under 
the action of G [3[T4]: 

,9o € Go : D(go)V = V . (31) 

Of course the embedding of Go within G depends in general on V ■ Let us show that the scalar fields 
ip a spanning the submanifold Gq/Hq, Hq being the maximal compact subgroup of Go, are flat directions 
of the potential and of the W- function, namely that neither V nor W, depend on ip a . Since we are 
interested in the part of the little group which has a free action on the moduli, we shall dchne Go modulo 
compact group-factors. For instance if the little group is SU(3) x SU(2, 1), we define G to be SU(2, 1) 
and thus H — U(2). For a summary of the orbits of regular extremal black holes in the various theories 
and of the corresponding moduli spaces Gq/Hq see Table 1. 

To prove that <p a are flat directions of both V and W, let us decompose the n scalar fields <j> r into 
the ip a scalars parametrizing the submanifold Go / Hq and scalars ip k , which can be chosen to transform 
linearly with respect to Hq. Let us stress at this point that the coordinates ip a , ip k will in general depend 
on the original ones <p r and on the electric and magnetic charges, namely: 

V a = <p a (<t> r ,p A ,qA) , <p k = <p k (<t> r ,p A ,q A ). (32) 

Let us choose, for convenience, a basis of coordinates in the moduli space such that the first n/ components 
of <fi r coincide with the (p a , the others being ip k , that is <p a = ip a , (f> k — ip . We can move along the (f> a 
direction through the action of isometries in Go- We shall consider infinitesimal isometries in Go whose 
effect is to shift the a-scalars only: 

g Q eG Q : cb r -> ( 5o *0) r = </> r + W Q , V -> V' = V + 5V = V, (33) 

where we have used the definition of Go, (|3ip . Let us now use eq.s ([2l)f and (|2"T|) to evaluate the 
corresponding infinitesimal variations of V and W: 

V{4> r ,V) = V((f> r + 64> r ,V + 5T) = V(4> k ,<p a +5(f> a ,V). 

W(4> r ,V) = W(<j) r + 5(j) r \V + SV) = W((f> k > (f> a + 6(j) a ,V). (34) 

We conclude that = ^pj- = 0, namely that (f> a are flat direction of both functions. Using eq. (|TTj) 
we see that the same property holds for the ADM mass: -^Madm = 0. Let us now give a general 
characterization of the VF-function in terms of the central and matter charges. We can write the coset 
representative ~L{(jf) of Ji sca i as the product of the Gq/Hq coset representative L (<^ a ) times a matrix 
Li(0 fe ) depending on the remaining scalars: 

L(<f> r ) = L(0 Q ,0 fe ) ^LoOnM^). (35) 

We can write L(0 r ) as a 2ny x 2ny matrix L(</> r ) M ^, where M is an index in the real symplectic 
representation, while N spans a complex basis in which the action of H is block-diagonal. We can obtain 
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L(^) t ') m J y from the coset representative in the real symplectic representation 'hsp{ ( t ) ' r ) M n using the Cayley 
matrix: 

L(^) = h Sv m^ wh ere ^=-LQ _\M. (36) 

The central and matter charges Zab, Zi of the theory can be arranged, together with their complex 
conjugates, in a (2ny)-vector Z^ defined as follows: 

( Z AB \ 
Zi 
Z AB 
V 'z 1 J 



Zm(^,V) = 



= -h{<p r ) T CV = -U((f> k ) T L {(j) a ) T CV , (37) 



Now we can use the property of Lo (</>") of being an element of Go in the symplectic representation, so 
that JLjCP = CL^V = CV and write: 

Z A (4> a ,^,V) = -U{4> k ) T CV = Z £l {Q,4>\V), (38) 

that is the central and matter charges do not depend on <f> a at all: 

w Zab = w Zl ="- m 

Let us now describe the effect of a generic transformation go in Go on the central charges. From the 
general properties of coset representatives we know that D(go) Lo(0 Q ) = Lo(<7o *4> a ) DQio), D(ho) being 
a compensator in Hq depending on go and <p a . Now, using the property that 4> k transform in a linear 
representation of Hq, we can describe the action of go on a generic point 4> as follows: 

D(g ) h(<j> r ) = D(g ) U{<p a )U{<p k ) = Mffo * 4> a ) D(h ) L x (^) D{h y l D(h ) = 

= U{g^<t> a )U{<t>' k )D{ho)^H9^<p r )D{ho), (40) 

where (j>' k is the transformed of <fi k by ho, and (go * <j> a , <fi' k ) define the transformed go * 4> r of cf> r by go- 
From (|40]) and the definition (|37|) we derive the following property: 



V.g eG : Z fy (go*<l> r ,V) = \D{ho)- T ]M N Z ^ ,V) = h Q ^ Z ^ ,V) , (41) 

where, to simplify notations we have denoted by ho * Z the vector [D(ho)^ T ]£ I N Z-&. Now consider the 
W function as a function of </> r and V through the central and matter charges Z^: 

W{f,V) = W[Z A (4> r ,r)]. (42) 

From the duality- invariance of W it follows that, for any go € Go we have 

W(<f, V) = W(g * <f, D(g ) V) = W(g *<f> r ,V). (43) 

Furthermore, using eq.s (|4lT) . (l42j) we find: 

W[Z{cjf,V)] = W(<p r ,V) = W(g * <f, V) = W[Z(g * (f> r , V)} = W[h * Z^(<j) r ,V)] . (44) 

The above equality holds for any go € Go and thus for any ho £ Hq- We conclude from this that W 
can be characterized, for a given orbit of solutions, as an Ho-invariant function of the central and matter 
charges. This is consistent with what was found in [3J [5] . Let us stress once more that we have started 
from a generic charge vector V , so that the definition of Go, and thus of Ho, is charge dependent. We could 
have started from a given Go inside G and worked out the representative To of the G-orbit having Go as 
manifest little group. In this case, by construction, the (4> a , 4> k ) parametrization is charge- independent. 
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M 



H 


orbit 


H n 




Ri 






I 


K 6(2) 
SU(2)xSU(6) 
^6(6) 
USp(8) 






8 


E 7(7) 
SU(8) 


III 


42 


1 + 27 


6 


SO* (12) 
U(6) 


I 

II 
III 


SU(4,2) 
S[U(4)xU(2)] 

SU*(6) 
USp(6) 


(4,2)_ 3 + c.r. 
14 


(6,l) +2 + (l,l)_ 4 +c.r. 
2 X 1 + 14 


5 


SU(5,1) 


I 


SU(2,1) 


2+3 + c.r. 


3 X 1_2 + c.r. 


4 


SL(2,M) w 
SO(2) A 
SO(6,n) 


I 

II 
III 


aO(4,n) 
SO(4)xSO(n) 

SO(6,n-2) 
SO(6)xSO(n-2) 

on(l 11 ~ SO(5,n-l) 
OW(l, Lj X SO (5)xSO(n-l) 


(4,n) 
(6,n-2) 

(l,l) + (5,n-l) 


2x[(l,l) + (l,n)] 

2 x [(1,1) + (6,1)] 

2 x (1,1) + (5,1)+ 
+(l,n-l) 


Q 
O 


SU(3,n) 


I 

II 


SU(2,n) 
S[U(2)xU(n)] 
SU(3,ra-l) 


( 2 > n )n + 2 + C - r - 


(l,n)_ 2 +c.r. 

fill -A- r r 

\ A > L )l-n + c.r. 


S[U(3)xU(n)| 


S|0(3)xU(n-l)] 




SU(l,n+l) 
U(n+1) 


I 

II 


SU(l,n) 
U(n) 


n n+ i + c.r. 


l_ n + c.r. 




SL(2,R) 
SO(2) A 

X SO(2)xSO(n + 2) 


I 

II 
III 


SO(2,n) 
o\J ) X o\J(n ) 

qon 11 SO(l,n+l) 

bU(i,i)x so(n+1) 


(2,n) 
l + (n + l) 


2x [(2,1) + (1,1)] 
3xl + (n + l) 




Sp(6) 

am 


I 

II 
III 


SU(2.1) 
U(2) 

SL(3,K) 
SO(3) 


2_3 + c.r. 
5 


1-4 + 3+2 + c.r. 
2x1 + 5 


2 


SU(3,3) 


I 

II 
III 


/SU(2,1)\ 2 


(2, l) 3i0 + (l,2) 0i3 + c.r. 
8 


(2,2) 1 _ 1 + (1,1)_ 22 + C .r. 
2x1 + 8 


S[U(3)xU(3)] 


V U(2) J 
SL(3,C) 
SU(3) 




SO* (12) 


I 

II 
III 


SU(4,2) 


(4,2)_ 3 + c.r. 
14 


(6, l) +2 + (l,l)_ 4 + c.r. 
2 X 1 + 14 




U(6) 


S|U(4)xU(2)] 
SU*(6) 
USp(6) 




B 7(-25) 


I 

II 
III 


E 6(-14) 


16+3 + c.r. 
26 


1+4 + 10_2 + c.r. 
2 X 1 + 26 




U(l)xE 6 


U(l)xSO(10) 

E 6(-26) 
F 4 



Table 1: Summary of regular, extremal black hole orbits in the various supergravities. The symbols I, 
II, III denote the ^-BPS, the non-BPS (I4 > 0) and the non-BPS (I4 < 0) orbits respectively. For those 
solutions with non-trivial moduli spaces ^ (j. e . Q Q non-compact), the representations Ro, Ri of H , 
see Appendix...., are given. The symbol "c.r." stands for conjugate representations. 
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A detailed analysis. Let us exploit now, for the BPS and non-BPS extremal, regular solutions, 
the symmetry properties of the W function discussed in the previous sections, to study general aspects 
of the evolution of the flat and non-flat directions. 

We start computing the Killing vectors associated with the Go-transformations and write the con- 
dition that W be Go-invariant in the form of differential equations. To this aim, we will first compute 
the general expression for the vielbein of ^i sca \, in the parametrization (|35p . Let us denote by {T A }, 
A = 1, . . . , dim(Go), the generators of Go- We can perform the Cartan decomposition of the Lie algebras 
and go generating G and Go, respectively, with respect to their maximal compact subalgebras f), f)o: 

= Si ® f) , go = J? © f)o • (45) 

Under the adjoint action of H , the space Si split into subspaces ^o an d Si\ transforming in the repre- 
sentations Ro, Ri of H . The non-compact generators {K^}, f = 1 . . . , n, of Si (the indices f, s label 
basis elements of the tangent space to the manifold) split into the generators {K a }, a, b — 
of Sio, belonging to the tangent space of the submanifold Gq/Hq, and the remaining n — rif generators 
{Kj} of Si\. The Lie algebra f)o of Hq is generated by {H„}, u = 1 ... , dim(i?o)- As far as the choice 
of the parametrization is concerned, for the BPS and non-BPS (J4 > 0) solutions, we choose the coset 
representative as follows: 

U4> r ) = U{^ a )U{4> k ) G e Ro ■ e Al , (46) 

that is L O (0 Q ) is an element of e^ = Gq/H and hi((j> k ) is an element of e^ 1 . This in particular implies 
that c/> a and 4> k transform in the representations Ro, Ri of Hq, respectively (see Table 1 for a list of these 
representations) . 

As for the non-BPS (I4 < 0) solutions, it is more convenient to adopt a parametrization of the coset 
which is different from (|46p . in which <f> k can be defined to transform linearly with respect to the whole Go- 
For this class of solutions, see below, we define cf) k to be parameters of a solvable algebra {sk \ = {s\, so}, 
generated by n — n/ — 1 nilpotent generators s\ and a Cartan generator sq- As we shall see, for the 
standard choice of the charges Vo, <p k consist in n — rif — 1 axions originating from the D = 5 vector 
fields and a dilaton describing the modulus of the internal radius in the D = 5 — > D = 4 dimensional 
reduction. 

We want to compute the components of the vielbein of ^ S cai in the basis (|4"6"|) . To start with, the 
Gq/Hq left-invariant 1-form reads: 

ftoOn = £0^0 = ^T A = #^n /(^)T4. (47) 

If we split the Go-generators T A into generators of Gq/Hq (K a ) and of Hq (H„), that is A — > (a, it), then 
£lop b ((f> a ) = V/3 h (4' a ) defines the vielbein of Gq/Hq. Moreover, let us introduce the left-invariant 1-form: 

Qi((j) k ) = L^dL-L = d(t) k V k k ((t) k )K k + connection , (48) 

where d<fi k Vk k will define the vielbein 1-forms along the directions K^, of the tangent space. H 

In terms of the above quantities, we can now compute the left-invariant 1-form of ^# SCQ ; in the basis 
(H|: 

n(cj) r ) = h- 1 dh = hr 1 n h 1 +h 1 - 1 dh 1 = n A (<j) a )hi 1 T A h 1 + n 1 {(j) k ) 

= d(j)> 3 nop A {(t) a )h 1 /{(l) k )Kf + d(l) k V k k K je + connection, (49) 

6 The reason why the left-invariant 1-form S7i in eq. (|48|) does not expand on the generators K a will 
be clarified in Appendix O 
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where we have written Li -1 T4 Li = Ci A r {(j) k ) K.f + compact generators. Similarly we will also write 
Lq" 1 T a Lo = CoA^i^") Tb- The non-vanishing components of the vielbein V r s of ^ sca i are now readily 
computed: 

v B h = fio/ffJAA/) , v^ = n p A (r)c 1A \4> k ) , VfeWfeV). (50) 

Note that for all regular extremal black holes, our choice of parametrization is such that the vielbein 
matrix has a vanishing off-diagonal block Vk a , see Appendix [Cl 

The non-vanishing blocks of the inverse vielbein V~ l f- r are: 

v- 1 /, v-y, v-y = -v- v v/ v-y , (51) 

where V _1 a ^, V~y are the inverses of the diagonal blocks V a b , Vk k , respectively. 

Consider now an infinitesimal Go-transformation go ~ 1 + e^T^, e A ~ 0, and write (go * <^) r ~ 
r + e A k r A {cj)). The Killing vectors k r A {(f)) are computed, in the parametrization (|46|). to be: 

fc^ = ro/lflAB^^Jv- 1 /. (52) 

The Go- invar iance of the W-function (W(go * 0,"P) = W(</>, P)) can now be expressed in the following 
way: 

dW - , n dW 

*A W = «. £ 1 ^)V-V^=0, (53) 

where we have used the property that £0.4 (ck) is non-singular. Using the expression of the vielbein, it 
will be useful to write the first-order flow-equations for the scalar fields in the following form: 



cjfV/ = e u V- lrs ^- I \ TJ lt .££. ( 54 ) 

' 4> r V r k = e u V- lkk §%b 



We shall illustrate the implications of the above formula in two relevant cases: The BPS solution and the 
non-BPS one with I4 < 0. 

The BPS black holes. For the sake of concreteness we shall consider the supersymmetric regular 
solutions (i-BPS) in the maximal theory M = 8, although our discussion is easily extended to non- 
maximal theories. In this case Go = E 6 ( +2 ) and H = SU(2) x SU(6) C SU(8) = H. With respect to the 
adjoint action of Ho, the coset space &, in the 70 of SU(8), splits into the subspaces &o = {K a } in the 
(2, 20) and £x = {K^,} in the (1, 15) © (1, 15) of H , according to the branching: 

70 -> (2, 20) ffi (1,15)8 (1,15). (55) 

The parametrization (|35|) amounts to the following choice of the coset representative: 

L = Lo(0 Q )L!(0 fe ) , L (f)e^» , Li(0*) e e* 1 . (56) 

Since the index a spans a SU(2)-doublet (a = (A, A), A = 1, 2, A = [mnp] — 1, . . . , 20, m, n,p = 1, . . . , 6), 
while k only SU(2)-singlets, being tfi k themselves SU(2)-singlets, the non vanishing components of the 
matrix Ci A r are: £i a b (<fi k ), £i u k (<fi k ). Consider now the implications of the Go-invariance of W , as 
expressed by eq. ([53jl . The Hq = SU(2) x SU(6)-invariance corresponds to the A — u component of the 
equation, and implies 
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The invariance of W under Go/i?o-t r a ns f° rm atioriS, on the other hand, implies, using (fSTj) and (|50|) : 

o = c la \4> k )v-\ r ^=c la h {4 ) k ) 
= Ci a \4> k ) 



v -i,dW lk dW 



dW 



dW 



c la \^)v-^— => — = 0, 



90° 



(58) 



where we have used eq. (|57|) and the property that the block £i a b (cf> k ) is non-singular. The above 
equation expresses the (/("-independence of W, which we had proven before in a different way. Finally, 
consider the evolution of the </> Q -scalars as described in ([54]) . From equation (|58| it follows that: 

(59) 



namely the flat directions cj) a are constant along the flow. This is consistent with the AT = 2 supersym- 
metry of the solution, since the variation of the fermions A mnp (the hyperinos in the J\f = 2 truncation, 
in the 20 of SU(6)) on the solution reads: 



SX mnp oc <j> a V a A ' mnp e A = , 
where, as usual, we have written a = (A, mnp). 



(60) 



As far as the non-BPS black holes with I4 > are concerned, the analysis is analogous to the BPS 
case illustrated above. 



Non-BPS black holes with J 4 < 0. In this case the little group G of the charge vector 
is the duality group of the five-dimensional parent theory (for the AT — 8 case Go = E 6 ( 6 )), so that 
the fiat directions (4> a ) spanning Gq/Hq are the five-dimensional scalar fields. We can use the solvable 
parametrization for ^ S cai by writing ^ sca i — exp(Solv), where Solv the solvable Lie algebra defined 
by the Iwasawa decomposition of G with respect to H. Let moreover SoIvq be the solvable Lie algebra 
generating the submanifold spanned by the flat directions: Go/ Ho = exp(Solvo)- 

In the solvable parametrization the moduli <p a are parameters of the generators s a of So/voEI. We 
can decompose the scalars (jf into </>" and 4> k by decomposing Solv with respect to Solvo- 

Solv = o(l, 1) 8 So/^o 8 R-2 , (61) 

where the o(l,l) generator sq is parametrized by the modulus ao of the radius of the fifth dimension 
and the abelian subalgebra R_2 = { s a} is parametrized by the axions er A originating from the five- 
dimensional vector fields and transforming according to the representation R of Go with 0(1, l)-grading 
+2 (in the maximal theory R = 27). The decomposition (|6ip originates from the general branching rule 
of G with respect to Go 

Adj(G) = 1 O Adj(G)®R_ 2 ®R +2 , (62) 

The non-flat directions (f> k therefore consist of ao and o~ A , which transform in a representation of Go . The 
following commutation relations hold: 

[Ta, s q ] = , [s , sa] = +2sa , [T4, sa] = -Taa S s s . (63) 

7 Note that, in contrast to the parametrization used for the other classes of black holes, neither <fi a 
nor the corresponding solvable generators s a transform, under the adjoint action of Ho, in a linear 
representation. 
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We shall write h 1 (cj) k ) = L(<j a ) e a ° s ° . 

Note that the coset parametrizations that we are using throughout this section, defined in eq. (|46[) . 
differ from the standard parametrization of ^# sca ;, which originates from the D = 5 — > D = 4 reduc- 
tion (like, for instance, the special coordinate parametrization of the special Kahler manifold in the 
Af = 2 theory). The standard parametrization corresponds indeed to the following choice of the coset 
representative: 

L(0 r ) = L(a A )e CTOSO Lo(0 Q ). (64) 

The prescription (|35[) . that we are using here, yields instead a different parametrization in which the order 
of the factors in the coset representative is different: h((f> r ) = Lo(0 Q ) L(cr A ) e aos " . The two parametriza- 
tions are related by a redefinition of the axions: 

a A = L-^on^, (65) 

where Loe A (0") is the matrix form of ho((j) a ) in the R representation: Lo(</> a ) _1 Lo(0 Q ) = Lo s A (4> a ) S A- 
The vielbein 1-forms d<jf V r r are defined, as usual, as the components of the left-invariant 1-form along 
the non compact generators Kg oc (s r + s]S). The non- vanishing components of the vielbein matrix V r r 
and of its inverse V _1 f r are readily computed to be: 

V a \4> a ) , V a A = ~e- 2 ° a V a \r)s b v k ^ , Va 4 = e- 2CT0 sj 1 , V 6 = 1 , 
V-V(0 Q ) , V-\ A = a aT1 A a s , V-Y^e 2 ^ , V^ = 1 , (66) 

where s a s A is the matrix form of the generator s a of SoIvq in the representation R. Consider now the 
Go-invariance condition on W, as expressed by eq. (|53l) and use the following property: 

Li^T^Li = T A - er 2aa T AT K = £/K f + compact generators. (67) 

After some algebra we find that the 7?o-invariance of W (component A = u of eq. (1531) ) implies : 

tf^f^ = 0, (68) 
while the invariance with respect to Gq/Hq (component A = a of the same equation) implies: 

dW 

w ' m 

that is W must be a-independent, as expected by other arguments. Let us note however that now the 
cj) a are evolving since: 

^ a V aa = e u SaS A a s ^^0, (70) 

since the right hand side represents the variation of W corresponding to an infinitesimal Go / H$ transfor- 
mation of cr A and W is invariant only with respect to iJo-transformations of er A (see equation (|68[) ). One 
can easily verify that the flow of the non-flat scalars (oo, er A ) is described by an a-independent dynamical 
system which has an equilibrium point for = — 0, at which, by virtue of ([7TJ|) . also a = 0. Indeed, 
using eq.s (|54[) and the explicit form of the vielbein matrix and of its inverse (|66p . we can substitute in 
the equations for <f> k the expression of d a and find for the non-flat directions the following equations: 

•A VI 4ff rAS , a A £ A T\ • u dW 

a = e (e 6 + s A s aF a a ) , a = e ^— . (71) 
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According to the above equations, the non-fiat directions a , Co evolve towards fixed values at the horizon 
which depend only on the quantized charges and solve the equilibrium conditions J^- = = 0- Only 
the flat directions can depend at the horizon on the values of the scalar fields at radial infinity, but this 
is not in contradiction with the attractor mechanism since the near horizon geometry only depends on 
the corresponding values of a A , <7o, through V or W. 

Let us finally give an example of the (0 Q , (/> fe )-parametrization in the STU model, in the case Ia{V) < 
0, and show that the central and matter charges do not depend on a. The STU model is a JV = 2 
supergravity with n = 6 real scalar fields (i.e. 3 complex ones {s,t,it} = {zi, Z2, Z3}) belonging to three 
vector multiplets. The number of vector fields is ny = 4. The scalar manifold has the following form: 

f SLj2JR)\ / SL^R)\ f SL(2J*)\ 

where each factor is parametrized by the complex scalars s = a[ — ie Vl , t — a' 2 — ie^ 2 , u = a 3 — ie V3 . 
The eight quantized charges transform in the (2,2,2) of the isometry group G — SL(2,R) 3 and in this 
representation the coset representative is the tensor product of the coset representatives of each factor in 
(1721) in the fundamental representation of SL(2,R): 

L(zi) - L 1 (z 1 )®L 2 (z 2 )®L 3 (z 3 ), (73) 

where each 2x2 matrix has the following form: 

In this case the a' A axions are nothing but a'^a^ag. The little group of the Ii('P) < orbit is Go = 
0(1, l) 2 . For generic charges, like for instance those corresponding to the DO, DA system (qo,p l ), the 
action of Go is rather involved and depends on the charges themselves. We can consider however, as 
representative of the same G-orbit, the charges corresponding to the DO — D6 system (p°, qg). In this 
case Go is parametrized by two combinations of the dilatons ipf. {<j) a } a =i,2 = {0 1 = (Vi — P2), <f> 2 = 
-^(<Pi + f2 — 2^3)}. According to the general prescription (|35l) . the part Lo of the coset representative 
depending on the flat directions (j) 1 , (j> 2 , should be the left factor of the product. This corresponds to 
bringing the diagonal dilatonic factor in (|74|) to the left and redefining the axion: 



where ai = a! i e~ Vi . The three complex scalar fields, in this new parametrization, read: 2, = e Vi (dj — i). 
The central and matter charges have the following form: 

vl,... 

[go + /e^ CT °(a 1 - i)(a 2 - i)(o 3 - i)] , (76) 



e 2 



±1 



(Jo 



Zi = ^-^[qo+p e^^( ai +i)(a 2 ~i)(a 3 -i)}, (77) 



VS., 



Z2 = — ^=-[ go +/e^ CTo (a 1 -i)(a 2 +i)(a 3 -i)], (78) 
2 V 2 



vl, 
e 2 



- - [go+P e V ^ CTO (ai"i)(«2-i)(a 3 +i)], (79) 

where <tq = ~^(<Pi + ¥?2 + <fi3)- We observe that none of the central and matter charges depend on the 
scalars {<fi a } — {(ft 1 , (j) 2 }, but only on the remaining scalar fields {<fi k }, k = 3, . . . , 6, defined as follows: 

{0 fc } = {0-0,04 saje-*"*}. (80) 
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The scalars </>" are then flat directions of any function of the central and matter charges, including V and 
W. 



3. Small black holes in the Af = 8 theory 

Af = 8 supergravity admits two orbits for "large" extremal black holes (one of which is 1/4-BPS and 
the other a non-BPS one) and three orbits for "small" extremal black holes (all of them BPS, preserving 
1/8, 1/4, and 1/2 supersymmetry respectively). 

Following the analysis of [14], the ADM mass for the three small orbits is given by the largest 
eigenvalue of the central charge matrix Zab- Its eigenvalues for 1/8 and 1/4 BPS solutions are given by 
the quartic and quadratic roots of the secular equation 

4 

I1( A - A <) = ( A > = A) > ( 81 ) 

i=l 

pi being the skew-eigenvalue of Zab when written in normal form. In particular we have: 

• For 1/8 BPS: Ai > A 2 > A 3 > A 4 

• For 1/4 BPS: \ x = A 2 > A 3 = A 4 

• For 1/2 BPS: A* = A V i = 1, • • • , 4. 

The five Af = 8 orbits preserve, respectively, the following symmetries: 



• large: 



1/8 BPS 
non-BPS 

1/8 BPS 
. small { 1/4 BPS 



1/2 BPS 



SU(2) x SU(6) 

US P (8) 
USp(2) x USp(6) 
SU(4) x US P (4) 
USp(8) 



The superpotential W, for all the BPS orbits, is given by the highest eigenvalue of the central charge 
matrix Zab, however one can also get small orbits from the large non-BPS orbit, in the limiting procedure 
Ii — > 0. Indeed, in this limit the non-BPS orbit becomes supersymmetric, the fraction of supersymmetry 
preserved depending on whether further constraints on I4 are imposed. For example, let us start with a 
non-BPS black hole with charges (p°,qo) turned on. It has I4 = — (po<7°) 2 and symmetry USp(8). The 
limit I4 — 0, obtained for poq° = 0, gives a 1/2 BPS black home which has the same USp(8) symmetry. 
For the most general W of a non-BPS configuration, as defined in [5], the I4 = limit just gives back eq. 
(EU) with A = W 2 . 



4. Small black holes in the Af = 4 theory 

The moduli space of Af = 4 supergravity is 

_ SO(l,2) SO(6,n) 
Jlscai - SQ(2) x S0(6) x SO(n) ' 

The presence of a non-simple U-duality group G — SO(l,2) x SO(6,n) makes the analysis more 
involved here than in the Af = 8 case, and requires to explicitly set the notations. 
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We consider static and spherically symmetric extremal black holes. In particular we address our 
study to the small orbits of the theory, corresponding to a vanishing value of the horizon area. They are 
identified by a vanishing quartic invariant of the U-duality group. 

The electric and magnetic charges span the representation (2, 6 + n) of the U-duality group SU(1, 1) x SO(6, n 

VI, a =1,2, A = l,---,6 + n, 

such that 

V\ = PA = ry AS p S , V\ = q A ■ (82) 

6 n 

Here t]at = diag(+, •••,+,—,•••,—) is the SO(6, n)-invariant metric. 

It is useful to introduce an SO(l, 2)-invariant tensor quadratic in the charges: 

Tat, = ^{pAqs - qAPs), 

and an SO(6, n)-invariant tensor, quadratic in the charges, in the adjoint representation of SO(l,2), 
obtained as follows: 

L a = \l a aP V a ■ 

where 7°^ = (1, — 03, ax) a and a 1 denote the Pauli matrices. The indices a in the 3 of SO(l,2) are 
lowered and raised with the metric r) a b = diag(+, — , — ). The explicit form of the 3-dimensional vector 
L a is given by: 

L a = ±( P 2 + q 2 ), Li = ±(p 2 -q 2 ), L 2 =p- q (83) 

where p 2 = p a pa, q 2 = q A q A , p- q = p A qA- 

The quartic U-invariant of the theory is given by 

h = e^V^sV^V^ =p 2 q 2 -(p- qf = L a L b Vab = 2T AE T AS . (84) 

The covariant tensors Tat arL d L a can be expressed in terms of derivatives of the quartic invariant, 
restricted to the adjoint representation of the two subgroups of the U-duality group as [2] : 



24 dV A dVf dV 2 



(85) 

(l,Adj(SO(6,n))) 



8(5 + n) 7 dV A dV 0A dp 2 



(86) 

(Adj(SO(l : 2)),l) 



For the subsequent analysis we shall use the expression of I4 in terms of the central and matter charges 
given in (I16p . This formula can be understood by noting that the vector 

L a = (5i,Re(5 2 ),Im(5 2 )) 

transforms in the 3 of SO(l, 2), being related to L a through the action of the coset representative of the 
S so(2)" > f ac t° r in the same representation. Therefore the 2 vectors L a and L a are in the same duality 
orbit (indeed they coincide in the origin of the s ^j^ factor). This is a relation between dressed and 
bare charges which will be very useful in the sequel. 



8 In this and in the following sections the indices r, a, a, A have different range and definition with 
respect to the previous sections. Their definition will be given as soon as they are introduced. 
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When J4 > 0, the sign of L° (and hence of Si ) has a U-duality invariant meaning. Indeed, I4 represents 
the norm of the vector L a . Using the terminology of the Lorentz group, a positive or null norm vector L a , 
being "time-like" or "light-like" respectively, has the sign of its time component invariant under SO(l, 2) 
transformations. Viceversa, if I4 < L a is space-like and the sign of its time-component has no invariant 
meaning. 

Exploiting the symmetries of the theory, the central charge matrix can always be reduced to the 
fzie \ 

normal form Z^g — > with skew eigenvalues zi,z 2 G K, while the matter charge- vector 



z 2 e , 

Zi = pie 1 ® 1 can always be reduced to a form where the first two components are p\ e 10 and p 2 , all the 
other entries being zero [141 [2] . In this normal frame Si and 5*2 take the simple form: 

51 = z 2 + z 2 - p 2 - p 2 (87) 

5 2 = 2z 1 z 2 -p 2 1 e 2ie -p 2 2 , (88) 

and the general expression for the quartic invariant in terms of dressed charges in the normal frame reads: 
h = {zi-z 2 f [{z 1 +z 2 f- 2 (p? +p|)] +2 jO 2( jO 2„2z 1 z 2 )(l-cos(2 0)) • (89) 

We shall however use a different normal form for the matter charges, in which central and matter charges 
appear in a more symmetric fashion. 

Z[\ J_ /l i \ fZt\ = J_ ( P ie ie +ip 2 \ = £iv f 
Z'J ~ V2 U -i) \Z*J V2 \pie i9 ~ip 2 J C \p 2 e- ip . 



(90) 



The phase j3 can be absorbed by a SO(2) C SO(rt) transformation which, in this new basis reads 
e^P \ 

. „ . Taking into account that (Z^) 2 + (Z 2 ) 2 = 2 e 2l,p pi p 2 , Si and S 2 take the following 
e 1 1 J 

form: 

Si = z 2 + z 2 ~pi-pl S 2 = 2(ziz 2 -pip 2 e 2iifi ). (91) 

If we start from Zj = pie ie ' , I = 1,2, and fix 9\ — 62 = 7r/2, using the SO(2) C SO(n) freedom, wc 
easily find the following relations pi — Pl ^ 2 and p 2 = Pl ^ 2 ■ The quartic invariant in this new normal 
form reads: 

h = (fi-%) 2 + (z 2 -z 2 ) 2 -2(fi + pi)(z 2 + z 2 ) + 8pip 2 ziz 2 co S (2v) (92) 

Notice the symmetry between z±, z 2 and pi, p 2 . Indeed this normal form can be easily obtained from the 
J\f = 8 central charges by identifying zi,z 2 , pi, p 2 with the moduli pi of the skew-eigenvalues of Zab- 

It is straightforward to verify that for generic pi, p 2 the stabilizer on the matter sector is SO(n — 2). 
This does not change if either one or the other of the two norms is put to zero. If, on the other hand, 
Pi = P2 the stabilizer is enhanced to SO(n — 1) and, finally, if the two matter charges are both zero we 
recover the full SO(n). The five parameters {zi, z 2 , pi, p 2 , ip} defining the normal form can be expressed 
in terms of the five iJ-invariant functions of the central and matter charges characterizing a generic 
configuration of scalar fields and charges at infinity, i.e. the if-orbit of the solution. The classification of 
small and large black holes, on the other hand, refers to the G-orbits of the quantized charges V . Each 
G-orbit will in general comprise infinitely many i7-orbits, defined by G-covariant conditions on the five 
invariant parameters. 

4-1- Large orbits 

Let us first consider the large black-hole solutions, with J4 ^ 0. We have three orbits, which can be 
characterized, all over the flow, in terms of the sign of I4 and Si: 



19 



a) 1/4-BPS orbit, preserving the symmetry S0(4) x SO(n). 
In this case I4 > 0, and Si > 0. 

(3) non-BPS orbit with Z AB ^ 0, preserving USp(4) x SO(n - 1). 
In this case I± < 0, and the sign of S\ has no restrictions. 

7) non-BPS orbit with Z^b = 0, preserving the symmetry SU(4) x SO(n — 2). 
In this case I4 > 0, and Si < 0. 

4.2. Small orbits 

Let us now consider the case of small orbits, corresponding to I4 = 0. They can be classified, in 
terms of the invariants introduced above, into 3 inequivalent classes |14j : 

A) T AS = 0. 

This class contains 3 different orbits, corresponding to L° being positive, negative or null: 

Al) 1/2-BPS orbit preserving the symmetry USp(4) x SO(n). 

In this case L°,Si > 0. 
A2) non-BPS orbit preserving the symmetry SU(4) x SO(n — 1). 

In this case L°,Si < 0. 
A3) 1/2-BPS orbit preserving the symmetry USp(4) x SO(n - 1). 

In this case L a = 0, Si = S 2 = 0. 

We note that this class of orbits has a simple realization in the heterotic basis, where the charges 
Pa = V\ = 0. In this basis the three orbits correspond to the norm of <ja being positive, negative, 
or null. 

B) T As ^0;L a = 0. 

This is a 1/4-BPS orbit preserving the symmetry SO(4) x SO(n — 2). 

C) 7W0;L a ^0. 

It contains two orbits: 

CI) 1/4-BPS orbit preserving the symmetry SO(4) x SO(n - 1). 

In this case L°,Si > 0. 
C2) non-BPS orbit preserving the symmetry USp(4) x SO(n — 2). 

In this case L°,Si < 0. 

The above classification was found by studying the near-horizon properties of the solutions. However, 
since the different orbits are characterized in terms of U-duality invariants, actually the same properties 
hold true all over the flow of the fields from space infinity to the horizon. 

4-3. The W function for small orbits 

In the study of the large orbits of extremal black holes, the dynamical flow can be completely 
characterized in terms of the fake-superpotential W , which enjoys the property of being a monotonic 
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function decreasing from space- infinity (at r = 0) and the horizon (at r — » — oo). In particular the ADM 
mass is defined in terms of W as: 

M ADM = lim U = lim Jw(^(t)) = W(0 O ) , (93) 

and it is bounded from below by the value taken at the horizon, which is fixed in terms of the charges 



M ADM = W(fo) > W{4>{-oo)) = \ -T-. (94) 

V 4 7T 

In this respect small orbits are problematic because in this case the mass is not bounded from below. This 
raises the problem of defining the fake superpotential for small orbits, where the attractor mechanism 
breaks down. On the other hand, since W is well defined for the large orbits, we may nevertheless try to 
define it from the large orbit cases by an appropriate limiting procedure obtained by choosing particular 
constraints equivalently on the bare or on the dressed charges, such that the horizon area collapses to 
zero: 

A H = 4tt -> 0. (95) 
This requires a careful analysis of the behavior under this limit of the invariants of the various orbits. 



4-4- Small versus large orbits 

In the following we give the constraints on the charges needed to obtain zero-horizon area for each 
of the orbits listed above. To this end let us describe the large orbits in terms of the five normal 
form parameters Zt, z%, pi, p2, (p. In particular we shall consider representatives of the G-orbits of the 
quantized charges and the corresponding values of the five parameters on some specific point 0o of the 
moduli space at infinity, keeping in mind that on a generic point the solution is characterized by five free 
parameters. 

a) Being the solution 1/4-BPS, in a generic point </>o at infinity the W function is given by the highest 
eigenvalue of the central charge matrix, say z%, expressed in terms of H- invariant functions of the 
central and matter charges: W = \zi\. As a representative of the orbit we can take p 2 = and 
denote p\ = p. In this case S 2 is real and thus: 

h = (Si~S 2 )(Si + S 2 ) = [p 2 -(z 1 +z 2 ) 2 ][p 2 -(z 1 -z 2 ) 2 }. (96) 

Requiring it to be positive together with Si leads to the condition p 2 < (z\ — z 2 ) 2 . As we shall 
see in the following, to obtain the A-type small orbits by setting some of the charges to zero, the 
above representative is not useful but a different one should be chosen, with p\ ~ p 2 . 

0) We can choose p~\ = p 2 = p and denote z\ = z 2 = z. In this case we have: 

St = 2(z 2 -p 2 ), 

S 2 = 2(z 2 -p 2 e 2i *) : 

h = -8p 2 z 2 (1 - cos(2< < 5)) < 0. (97) 

For this orbit the general expression of W in terms of _ff-invariant quantities is not known; 

7) The W function for the (7) large orbit was computed in [2]. There it was shown that, fixing the 
relative phase of the two matter charges to it/2, W was simply expressed as W = (P1+P2) = Pi, 
no constraints being required on the invariant quantities zt, z 2 , pj, tp characterizing the solution 
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at infinity. Just as in the (a) case, a representative of this G-orbit can be chosen by setting p 2 = 0, 
pi being denoted by p. The quartic invariant is given by eq. (|96l) and the condition I4 > 0, Si < 
requires taking p 2 > (zi + z 2 ) 2 ■ 

This representative, however, is not useful to retrieve the A-type small orbits by simple constraints 
on the remaining charges and a different representative should be chosen (having p\ = p 2 ). 

The large orbit (J3) has no constraints on Si , so it can generate all the small orbits by imposing appropriate 
relations among the charges. However, this orbit is the only one for which we do not have a complete 
knowledge of W, except for those ii-orbits with fixed <p. We will then be interested in deriving the small 
orbits from simple conditions on the charges describing the (a) and (7) large ones, for which W is known. 
The W functions of the small orbits will then be obtained from those describing the (a)- and (7)-solutions 
by imposing appropriate constraints on the quantized charges. 

The A-type orbits. As pointed out earlier, the condition Tae = can be solved by setting p\ = 0. 
Under this condition, as shown in Appendix |DJ central and matter charges satisfy a reality condition, 
which, on the two normal forms, implies: 

{z\ = z 2 = z , 9 = 0} {zt = z 2 = z , pi= p 2 = P , <P = 0} . (98) 

We can easily see from ([92]) that the above conditions imply I4 = 0. In particular, the condition ip — > 
implies 

U -> [{zi - z 2 f - (pi - p 2 ) 2 } [(zi + z 2 f - (pi + p 2 ) 2 } (99) 

Al) Since this orbit has S\ > 0, it can be derived from the large orbits (a) and (J3). 
The orbit (a) in the normal form has: 

z 2 1+ z 2 2 > pj+p 2 2 . (100) 

To obtain the G-orbit (Al) from the (a)-orbit it suffices to set Z\ — > z 2 , pi = p 2 = consistently 
with (f9"5)) . Then from I4 — > 0. The compact symmetry of the charges is enhanced to Hq = 
USp(4) x SO(n) and the preserved amount of supersymmetry is doubled to 1/2-BPS. 

The W function is given, in a generic point on the moduli space at infinity, by the highest eigenvalue 
of the central charge matrix: W = z\ ; 

A2) It has Si < 0, so that it can be obtained from orbits (/?) and (7). This is consistent with the fact 
that it is a non-BPS orbit, which cannot be obtained from a BPS one by imposing relations on 
the charges. The orbit (7) in the normal form has: 

z 2 + z 2 < p 2 +p 2 . (101) 

To obtain a representative of the G-orbit (A2) from (7) we can start from a representative of the 
latter with z\ = z 2 = and impose condition (|98[) which requires setting pi = p 2 . 

The compact stabilizer of the charges is enlarged, with respect to the original orbit (7), to H = 
SU(4) x SO(n- 1). 

The W function is given, in a generic point on the moduli space at infinity, by the highest of the 
two matter charge moduli pi: W = p±: 
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A3) This ^4-orbit has Si = 0. It can be obtained from all the large orbits (a), (J3) and (7). In particular, 
it is found from orbits (a) and (7), by imposing (|98|) . so that in particular Si reduces to 

Si = 2(z 2 - p 2 ) (102) 

supplemented by the condition 

z 2 = p 2 (103) 

In this limit the compact stabilizer of quantized charges reduces to Hq — USp(4) x SO(n — 1). 
Being the solution 1/2-BPS, the W function is given, in a generic point on the moduli space at 
infinity, by the highest eigenvalue of the central charge matrix: W = z\ . Note that eq. (|103[) is 
a condition between central and matter charges, and it is a necessary condition, when the orbit 
is obtained from the (7) large orbit, for the enhancement of supersymmetry from non-BPS to 
1/2-BPS, since for the (7) orbit W is given by the highest eigenvalue of the matter charges, while 
in the BPS cases it is given by the highest eigenvalue of the central charges. 



The B and C type orbits. These orbits are characterized by having T AS ^ 0. This implies that 
some of the conditions should be relaxed. In particular we could have the following possibilities: 



B) 
CI) 
C2) 



Zi ^ z 2 
z\ ^ z 2 
Zl = z 2 



P\=P2, 

pi + Pi- 



(104) 



Let examine these three cases in some detail. 



B) The whole vector L\ a is zero in this case, which implies S\ = S 2 = 0. We can start from a 
representative of either the (a) or the (7)-orbit with z 2 = = p 2 . In the former case z\ > p\ while 
in the latter zf < p 2 , while Si and I4 read: 

Si = z\ - p\ , S 2 =0, 

h = S 2 = (z 2 -p 2 ) 2 . (105) 

To obtain a representative of the £?-orbit we need to set: 

z\ = p\ (106) 

The compact little group of this orbit is Hq = SO(4) x SO(n — 2), and, being the solution 1/4-BPS, 
the W function on a generic point of the moduli space is given by the highest eigenvalue of the 
central charge matrix, say W = z%. Similarly to the (A3) case, we note that the condition (jl06[) 
between central and matter charges is necessary, when the orbit is obtained from the (7) large 
orbit, for the enhancement of supersymmetry from non-BPS to 1/4-BPS. 

CI) It has Si > so it can be obtained from orbits (a), (/?). We can start from a representative of (a) 
with ip = and pi — p 2 = p (but zi > z 2 ). In this case 

h = (zi-z 2 ) 2 [( Zl +z 2 ) 2 -4p 2 ]>0. (107) 

If we further impose 4p 2 = (zi + z 2 ) 2 , we find I4 = and Si = h (zi — z 2 ) 2 > 0. The compact 
little group of this G-orbit is H Q = SO(4) x SO(n - 1) and, being it 1/4-BPS, the W function is 
given, in a generic point on the moduli space at infinity, by the highest eigenvalue of the central 
charge matrix: W — zi . 
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C2) It has Si < so it can be obtained from the orbits (f3), (7). In the latter case we can start from 
a representative with z\ = Z2 = z, p2 = and p = pi > 4z 2 . In this case Si and I4 read: 

Si = 2z 2 -/5 2 <0, h = P 2 (P 2 -4z 2 ) > 0. (108) 

A representative of C2 is obtained by setting p 2 = 4z 2 . 

On a generic point in the moduli space, the W function is given by the highest of the two matter 
charge moduli pi, say W = p\. 

5. Small Black Holes in M = 2 Magic Models 

We are going to show that, as in the Af = 4 case, the M = 2 ADM mass for small black holes either 
is supersymmetric (and then given by |Z(0o)|, in terms of the asymptotic value at r — > of the central 
charge) or can be obtained by a non-BPS black hole with Z = 0, and it has therefore a known expression 
in terms of radicals [XSJ E] . An exceptional case in this respect is the t 3 -model which has no Z — orbit. 
In this case, as we shall show below, the W for the small orbit can be obtained either from that of the 
BPS I4 > orbit or from the VF-function of the non-BPS I4 < orbit, which is known [3]. 

5.1. D = 6 Uplift of Magic Models 

In D = 6 there are three types of black holes, corresponding to the sign of the quadratic form 

X 2 = X A VAB X B , (109) 

written in terms of the scalars in the tensor multiplets X A 6 SO(q+ 1, l)/SO(q+ 1) (q = 1, 2, 4, 8 for real, 
complex, quaternionic o octonionic models respectively). In particular, for X 2 ^ we have large black 
holes, associated to two-charge solutions, BPS for X 2 > 0, generated by two quantized charges of the 
same sign, and non-BPS for X 2 < 0, corresponding to quantized charges of opposite sign. For X 2 = 
we have instead a small black hole, which is a BPS one-charge solution. 

The compact symmetry of the three different orbits are 

• SO(q + 1) for the BPS large black hole 

• SO(q) for the non-BPS large black hole 

• SO(q) for the BPS small black hole 

We are going to consider in the following the octonionic case (q = 8). 

5.2. D = b Uplift 

There are three orbits of small black holes at D = 5, corresponding to the vanishing of the cubic 
invariant -Z3 = 0. Two of them are two-charge configurations, BPS if the 2 charges have the same sign, 
non-BPS otherwise, while the third is a BPS one-charge configuration, corresponding to -Z3 = 19/3 = 0. 
These three orbits are the trivial dimensional reduction of the six-dimensional orbits. They preserve the 
compact symmetries SO(9), SO(8) and SO(9) respectively. 

Moreover, there are two large orbits corresponding to three-charge configurations (with ^3 ^ 0), a 
BPS one for charges of the same sign, and a non-BPS one, when one of the charges has opposite sign 
with respect to the other two. The compact symmetries preserved are F4 and SO (9) respectively. 
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5.3. D = 4 Analysis 

At D — 4 there are five small orbits. Three of them originate from dimensional reduction of the small 
D = 5 orbits, and have compact symmetry 0(10), 0(9) x 0(2) (critical, two charges) and F4 (double 
critical, one charge) ; two of them are three-charge light-like orbits coming by direct dimensional reduction 
of the large D = 5 orbits and have compact symmetries F4 and SO (9). 



5.4. ADM Mass for Small Black Holes at D = 4. 

For the three BPS orbits it is given by the asymptotic value, at space infinity (r — > 0) of the norm of 
the central charge Z, at I4 — 0. For the two non-BPS orbits, both can be obtained by a large non-BPS 
Z = orbit (with charges (+,+,—,—)) by setting one of the charges to zero. Then the ADM mass is 
given by the asymptotic limit of W , which is obtained from the one of the corresponding large orbit, 
which is a known radical function of the H-invariants (see eq.s (5.11) and (5.12) of [5] and eq. (3.4) of 
0)- 

5. 5. The W Function of the Small Black Holes in the t 3 -Model 

As anticipated, the t 3 -model deserves a spacial treatment in our analysis since it does not have the 
non-BPS orbit with Z = 0. However, in this case, both the W functions for the I4 > (BPS) and the 
I4 < (non-BPS) orbits are known. The former is W = \[T\ = \Z\, while the latter is given in eq. (3.15) 
of [1] and reads: 



W 2 = 



i 1 + i 2 




(110) 



It is straightforward to verify that the above expression has a finite Zj — > limit which yields W — \fi\ = 
\Z\. In the J4 — > limit therefore, both the W- functions associated with the two large black hole orbits 
coincide with the modulus of the central charge, as expected since small black holes in this model are 
BPS. 



A Proof of eq. ([15]) 

Equation (| 13[> . as shown in 8|, is a particular form of the general solution to the Hamilton- Jacobi 
equation. In what follows we shall tailor the formal proof given in J3J to the class of extremal solutions 
we are considering, without making use of the Hamilton- Jacobi formalism. 

Consider the extremal solutions U(t; <j>o,Uo) and 4>(r;^>o,Uo), for a given charge vector V, within 
the interval < r < to, where now £/o, <fio denote the values of the fields computed at tq: Uq = 
U{t ; 0o, Uq), 0o = 4>( T o'i 4>o,Uq). The values of the fields at r*, for our family of solutions, is completely 
fixed in terms of (Uq, <f>o) and V . Let us perform an infinitesimal variation of the boundary conditions: 
Uq ^ U~o + SUq and <pQ = 4>o + 5<fio. This will determine a new solution within the same class: 

U(t; 0o + 8<h,Uo + 5Uo) = U(r; O) U ) + 8U(r) , 

<l>(T;<h+S<h,Uo + 8U a ) - 0(r; O , U ) + 8</>(t) . (Ill) 
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Now we write a seemingly more general ansatz for W than the one in eq. (fl"3 



e Uo W(^V) = e u ' W((j>*,V) + J e 2U ^ Uo ^ V(4>(t; U , <fo),V) dr . (112) 

As we shall see, the result of this integral does not depend on the choice of To. For the sake of simplicity 
we shall suppress the dependence on r and on the boundary values of the fields in the integrand. Since 
the integral is computed along solutions, we can use the Hamiltonian constraint (|6|) to rewrite W as 
follows: 



W(ct>o,V) 



e u * W(<f>*,V) + 



1 



1 



TO 



1 



V((p, V) + U z + - G r 



dT 



W(^,V) + - / C ef f(U,4>,U,4>)dr. 



Now perform the variation fj 1 1 1[) . integrate by parts and use the equations of motion: 



SU e u ° W(4>q,'P) + e u ° drW^o.V) > 



5(e u > W(fa,V)) + 



d_ r ±d_ r 

dU eff drdil ff 



SU 



1 



+ (USU + -G rs (p s 5<P) 



d_ r ±JLr 

W eff drrfidf eff 



1 



= 6(e u * ,V))+ (USU + -G rs ^ HI 



(113) 



where we have used the short-hand notation d r W = §^r. We can choose = — oo, so that all terms 
computed at in the above equation vanish. Equating the variations at To on both sides we find: 

U(t q ) = e Ua W{(j> 0i V) , <p s (To) = 2e Uo G rs (<f>o)d r W(<f>o,V). (114) 

Being tq generic, we find that W defines the first order equations ([7]) for the fields and thus it is a solution 
to the Hamilton- Jacobi equation. 

Note however that W, as defined in (|112[) . may in principle depend on the chosen value of To, that 
is W = W(Uq, <fio, To, V). Let us show that this is not the case, namely that W(Uq, 4>o, t o + St,V) — 
W(Uq-, 00; 7"0i for a generic St. To do this we vary To — > To + St, keeping the boundary values of the 
fields fixed. This requires to change the solution on which the integral is computed from [/(t),0(t) to 
U'(t), </>'(t) such that: 

U'(t + 5t) = U (t ) = U (t + St) -U (t ) St, 

4>'(to + 5t) = 0(t o ) = <b(T + St) ~ <j>(To) St . (115) 

and thus amounts to performing, along the flow, the transformation U St, <f) 4> — 4> St, 

besides changing the domain of integration, St being chosen along the flow so that St* = 0. After some 
straightforward calculations we find: 

e u ° (W(U , 0o, T + St, V) - W(U , fa, t q , V)) = H eff \ To St = 0, (116) 

in virtue of the Hamiltonian constraint. Since the function W of the moduli space, as defined by (|112p . 
does not depend on the choice of To, we can choose To = 0, where Uq — and then find ([T3]) . 



B Stability and Asymptotic Stability in the Sense of Liapunov 

Let us briefly recall the notion of stability (in the sense of Liapunov) and of attractiveness of an 
equilibrium point. Given an autonomous dynamical system: 

r = r (H7) 
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an equilibrium point </>* (f r (4>*) — 0), is attractive (or an attractor), for r — > — oo, if there exist a 

neighborhood X<^ of (f>*, such that all trajectories <j) r {r, (f>o) originating at r = in (f>o G evolve 
towards </>* as t — > — oo: 

lim r (r,0 o ) = 0: , V0 o el 0t . (118) 

T— >■ — OO 

An equilibrium point (not necessarily attractive) is stable (in the sense of Liapunov) if, for any e > 0, 
there exist a ball Bs{4>*) of radius S > centered in <^>*, such that: 

V0oGB ff (^). Vr<0 : 0(r,0 o ) e6 £ (0*), (119) 

that is, provided we take the starting point ^0 sufficiently close to <j)*, the entire solution will stay, for all 
t < 0, in any given, whatever small, neighborhood of 0*. Finally an equilibrium point is asymptotically 
stable (in the sense of Liapunov) if it is attractive and stable. 

Liapunov' s Theorem: If there exist a function v(4>) which is positive definite in a neighborhood of 0* 
(that is positive in a neighborhood of 0* and u(</>*) = 0) and such that also the derivative of v along the 
solution, in the same neighborhood, is positive definito ^ = (f> r d r v > 0, then (f>* is an asymptotically 
stable equilibrium point or, equivalently, a stable attractor. 

For large extremal black holes such function is v(<f>) = W((f) — W(4>*) = W(4>) — \f\Ii\- 

C Properties of the Vielbein on ^# scct z 

Let us briefly motivate why, for all regular extremal black holes, our choice of parametrization is such 
that the vielbein matrix has a vanishing off-diagonal block Vk a . 

The reason is purely group theoretical. As far as the BPS and non-BPS (I4 > 0) orbits are concerned, 
taking into account that <f> k belong to Ri and the index a label the Ro representation, d(j> k Vk a {4> k ) can 
be different from zero only if Ro is contained in the tensor product of a number of Ri representations. As 
the reader can ascertain from Table 1, this is never the case. For example in the case of regular BPS black 
holes, for N > 2, R is a doublet with respect to an SU(2) subgroup of H , while Ri is a singlet with 
respect to the same group. If we think of the JV = 2 truncation of the original theory of which the same 
black hole is a 1/2-BPS solution, this SU(2) group is the quaternionic structure of a quaternionic Kahlcr 
submanifold of the scalar manifold spanned by the scalars 4> a . On the other hand, as far as the non-BPS 
solutions with I4 < are concerned, the above argument does not apply in the coset parametrization 
(|46)l . but choosing instead the solvable parametrization one finds d(j) k Vk a (4> k ) = since Li -1 dL\ belongs 
to the same solvable algebra spanned by <p k , which is orthogonal to $0. 

D An Explicit Parametrization of the Af = 4 Scalar Manifold 

The quantized charges V = {V M ) of the Af = 4 model transform under G = SL(2,K) x SO(6, n) 
in the representation (2,6 + n), which can be labeled by the couple of indices M = (a, A), a = 1,2, 
A = 1, . . . , 6 + ?i. In this basis the symplectic invariant matrix reads: 

CmN = C a A : — — £ Q( 3 7iAE , (120) 

9 Here we require positive definiteness because our critical point is located at r — > —00 and not at +00 
as in standard textbooks. 
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Identifying the magnetic and electric charges with the components: V 1 A = p A , V 2A = q A = ?7 AS q-%, cor- 
responds to choosing the symplectic frame originating from the toroidal compactification of the heterotic 
string, in which SO(6,n) acts block-diagonally. We shall denote by L2 = £2 A s the coset representative 

° f SO(6)xSO(n) aild ^ Ll = tnat ° f S S0 2 (2^ ' denned as: 



la" P = (Xl el)' < m > 

where S = A — i is the complex scalar parametrizing the lower- half plane ^q 2 ^ • 

Using the definition (|37j) of the central and matter charges, we can write Z^, as Z a A, where a = 1, 2 
labels the complex basis in which SO(2) acts diagonally and which is obtained from the real one through 

the action of the Cayley matrix A = {A a ^) = .V In this notation Z~ 1K = (Z r , Zj), where 

Z r = 5 T r AB Zabi r = 1, . - . , 6, are the central charges and Zj, I = 1, . . . , n, are the matter ones. We 
find: 

ZaA — Aa 13 Ll" p L2 1 KiaSVY , (122) 

Using the explicit parametrization (I121[) we find: 

Z iA = (Z r ,Z I ) = 1 L 2 r A (q r + Sp r ) . (123) 

y/-2im[S) 

We see that (Z r , Z{) are real iff p A = 0. Reality of the Z r implies the following reality condition on the 
Zab- 

[ CD 

Zab — 2 € abcd Z . (124) 

If Zab is in the normal form, the above condition implies z\ = z%. Similarly, if the non vanishing matter 
charges ^7=1,2 = Pi: characterizing the normal form, are real, from (|90j) we find that Z[ — (pi+i P2) = 
pi e 1 ^, Z' 2 = Z[ = f>2 e -1 ^, and thus that p\ — p2- 
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